Quantification and scaling of multipartite entanglement in continuous variable systems 
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We present a theoretical method to determine the multipartite entanglement between different partitions of 
multimode, fully or partially symmetric Gaussian states of continuous variable systems. For such states, we 
determine the exact expression of the logarithmic negativity and show that it coincides with that of equivalent 
two-mode Gaussian states. Exploiting this reduction, we demonstrate the scaling of the multipartite entangle- 
ment with the number of modes and its reliable experimental estimate by direct measurements of the global and 
local purities. 
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The full understanding of the structure of multipartite quan- 
tum entanglement is a major scope in quantum information 
theory that is yet to be achieved. At the experimental level, it 
would be crucial to devise effective strategies to conveniently 
distribute the entanglement between different parties, depend- 
ing on the needs of the addressed information protocol. Con- 
cerning the theory, the conditions of separability for generic 
bipartitions of Gaussian states of continuous variable (CV) 
systems have been derived and analysed 111 yl |3|] . However, 
the quantification and scaling of entanglement for arbitrary 
states of multipartite systems remains in general a formidable 
task ^. In this work, we present a theoretical scheme to 
exactly determine the multipartite entanglement of generic 
Gaussian symmetric states (pure or mixed) of CV systems. 

We consider a CV system consisting of N canonical 
bosonic modes, associated to an infinite-dimensional Hilbert 
space and described by the vector X of the field quadrature 
operators. Quantum states of paramount importance in CV 
systems are the so-called Gaussian states, i.e. states fully char- 
acterized by first and second moments of the canonical op- 
erators. When addressing physical properties invariant un- 
der local unitary transformations, one can neglect first mo- 
ments and completely characterize Gaussian states by the 
2N X 2N real covariance matrix (CM) <j, whose entries are 
a,j = l/2{{X^,Xj}) - {X,){Xj). The CM a must fulfill 
the uncertainty relation cr + ift > 0, with the symplectic 
form V, = ®2^iUJ and uj = %-i — Sij+i, i,j = 1, 2. Sym- 
plectic operations (i.e. belonging to the group Sp(2N.'R) = 
{S e SL{2N, R) : S'^nS = fl}) acting by congruence 
on CMs in phase space, amount to unitary operations on 
density matrices in Hilbert space. In phase space, any iV- 
mode Gaussian state can be written as cr = S'^i^S, with 
i^ = diag{ni,ni, . . . nATjUAf}. The set E — {n.;} consti- 
tutes the symplectic spectrum of a and its elements must ful- 
fill the conditions rii > 1, ensuring positivity of the density 
matrix g associated to cr. The symplectic eigenvalues n^ can 
be computed as the eigenvalues of the matrix |iJ7cr|. They are 
determined by N symplectic invariants associated to the char- 
acteristic polynomial of such a matrix: two global invariants 
which will be useful are the determinant Dct c = Jli "^f ^^^ 
the seralian A = J^i iT'h which is the sum of the determinants 



of all the 2x2 submatrices of a related to each mode. 

The degree of mixedness of a quantum state g is charac- 
terized by its purity fi = Ti g^. For a Gaussian state with 
CM cr one has simply /i = l/V Deter. As for the entangle- 
ment, we recall that positivity of the partially transposed state 
g, obtained by transposing the reduced state of only one of 
the subsystems, is a necessary and sufficient condition (PPT 
criterion) of separability for {N + l)-mode Gaussian states of 
1 X A^-mode partitions jsllal- In phase space, partial transpo- 
sition amounts to a mirror reflection of one quadrature asso- 
ciated to the single-mode partition. If {n^} is the symplectic 
spectrum of the partially transposed CM cr, then a {N + 1)- 
mode Gaussian state with CM cr is separable if and only if 
iii > l\/i. A proper measure of C V entanglement is the loga- 
rithmic negativity Ej^/ IJ, which is readily computed in terms 
of the symplectic spectrum n^ of cr as Ej^ = — Yli-h <i In Jij. 
Such a measure quantifies the extent to which the PPT condi- 
tion ni > 1 is violated. 

Let us first consider the 2 A^ x 2 A CM cr^jv of a fully sym- 
metric A-mode Gaussian state (i.e. a state invariant under the 
exchange of any two modes) 
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where /3 and £ are 2 x 2 submatrices. Due to the symme- 
try of such a state, /3 and e can be put by means of lo- 
cal (single-mode) symplectic operations in the form (3 = 
diag{6, b}, e = diagjei, 62}. The symplectic spectrum 
Sniv of craN has then the structure (see the Appendix) 



{i/_, ...,iy-, iy+(N)}, 



(2) 



N-l 



vl = (6-ei)(6-e2), 
,^l^^, = (6+(A-l)ei)(fe+(A-l)e2). 

The (A — 1) -degenerate eigenvalue ;/_ is independent of 
A, while i^^(jv) can be expressed as a function of the pu- 
rity fip = (Det/3)~^/^ of the single-mode reduced state 



and of the symplectic spectrum of the two-mode block 0-^2 , 

E^2 = {jy_, v+ = i/+(2)} 

The global purity of the fully symmetric state is 



where i^_ is the lowest symplectic eigenvalue of the reduced 
two-mode state cr^2. The eigenvalues nip can be evaluated 
observing that Eqs. J4I7I9> impose the identity Acr = Ac,^ + 
{N — l)v^ + {v_ ^^fipN)~^ which can be used to obtain 
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and, through Eq. (|3}, can be directly linked to the one- and 
two-mode parameters. In particular, the symplectic eigen- 
values Vzf: are determined in terms of the two 5*^(4 r) in- 



Since partial transposition leaves the A^-mode symmetric 
block (TpN unchanged, the symplectic eigenvalues of a are 



variants /i^2 



and A«2 by the relation f^: 2v'i = A«2 =f ^8^*" "^ *^ f°™ ^ = {"'J' ^{v-, ■■■,^- 
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Next, we consider the [N + l)-mode Gaussian states con- 
stituted by generic single-mode states with CM a and fully 
symmetric A^-mode states with CM cr^w of the form Q. The 
mode with CM a is then coupled with all other N modes by 
the same 2x2 real matrix 7. The CM cr of such (7V+ l)-mode 
states reads 



r C/3-W 



(5) 



We will now show that the properties of mixedness and entan- 
glement of these states are determined by a suitable, limited 
set of global and local invariants under symplectic (unitary) 
operations. Let us introduce the purity /ia = (Dot a)^^/^ of 
the single-mode party, the global purity /ig. = (Dctcr)^^/^ 
of the state Q, and the global Sp(2N+2.'R) invariant A^- = 
Ac-y + A^jv — X^i 



tic spectrum E = {ni 



where the n^'s constitute the symplec- 
. . . , TiTv+i} of the CM cr, and 



Aq.^ = Det a + 27VDet 7 , (6) 

A^iv = iV(Det^ + (iV-l)Dct£) . (7) 

We are now in the position to characterize and quantify the 
bipartite entanglement between the single mode a. and the iV- 
mode block CTpN , the multipartite entanglement between all 
the A^ + 1 modes, and to provide an operational scheme for 
their experimental determination in terms of measurements of 
the global and local purities. To proceed, we must evaluate the 
logarithmic negativity by determining the partially transposed 
CM 6-, with respect to the partition a\j3^ , which is obtained 
by flipping the sign of Dct 7. Mixedness and entanglement 
are encoded respectively in the symplectic spectrum of cr, and 
of a. It is worth noting that, of the previously introduced 
parameters, only Aq,-,, is affected by the operation of partial 



transposition: Ac 



, with 



A„^ = Det a - 2iVDct 7 = - A„^ + 2/^^ 



(8) 



The symplectic spectrum E = {?ii} (i = 1, . 
CM (T Eq. Q is of the form (see Appendix) 



,iV+l)ofthe 



(9) 



defined as in Eq. MQ\ . but with A^-y replaced by A^^ from 

Eq. (|8j. The logarithmic negativity E'^ , quantifying the bi- 
partite entanglement between a and crpN, is determined only 
by those symplectic eigenvalues of a which satisfy ni < 1. 
Since v^ >1 (because it belongs to the symplectic spectrum 
of cr), the entanglement is determined only by the eigenvalues 
n^ . On the other hand, the eigenvalues n^ of Eq. dlOt can be 
interpreted as the symplectic spectrum of an equivalent two- 
mode state of CM er*^"? with global purity /i'^'' and serahan A^* 
given by 
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The corresponding A'^' associated to the partially transposed 
CM a"'^ reads then K<"i = -A'^? + 2/fil + 2/(j/^"V/3")^- 
By comparison with the expression A = —A + 2/iJ.f + 2//i2, 
holding for a generic two-mode state with local purities /ii 
and /i2 1 8], one determines the local purities of the equivalent 
two-mode state cr'^'^: 



eq 

/^l — Ma 



eq N- 

M2 = V- 



flf^l. 



(12) 



The two global invariants Eq. (Ill> and the two local invari- 
ants Eq. ( I12> determine uniquely the entanglement of the two- 
mode Gaussian state with CM cr"^"^. In particular, it is immedi- 
ate to see that the symplectic eigenvalues of the partially trans- 
posed CM a-'^'^ coincide with n^, so that we obtain the crucial 
result that the logarithmic negativity of the equivalent two- 
mode state coincides with the logarithmic negativity E^ 
of the {N + l)-mode state. Explicitely, one has: 



E 






maxjO, — logri-} , 



(13) 



with 2 711= A""} ~ JA^'"! -4:/^^<"!. Indeed, only the 
smallest symplectic eigenvalue n_ enters in the determination 
of the multimode entanglement, since n+ > 1 for two-mode 
states |8|. 

The 1 X A^ entanglement is completely quantified by mea- 
suring the two local purities jia and /i^n, the global purity 
/icr, the symplectic eigenvalue j/_, and Dct 7 (which together 
with /ic, determines Aa-,)- The experimental determination of 
these five quantities requires the full homodyne reconstruction 
of the (Af+l)-modeCMEq. (|5}. On the other hand, the study 



of the entanglement of two-mode Gaussian states has shown 
that a rehable quantitative estimate of the logarithmic negativ- 
ity, yielding exact (and very narrow) lower and upper bounds 
on the entanglement, can be obtained by simply measuring 
the global and local purities of the state 1 8]. In the present in- 
stance, this fact implies that a reliable estimate of the 1 x iV en- 
tanglement does not require the knowledge of the correlation 
matrix 7, while the remaining four quantities (the three puri- 
ties and the eigenvalue i/_) can be measured even without ho- 
modyning by direct single-photon detections |9|. Moreover, 
knowledge of these few quantities is also sufficient to deter- 
mine the multimode, multipartite entanglement of the state <j. 
In fact, the fully symmetric A^-mode block cr^w can be again 
regarded as a state describing a mode with CM (3 coupled with 
a fully symmetric [N — l)-mode block cr^jv-i, and thus the 
1 X {N — 1) entanglement within crpN can again be computed 
by constructing the corresponding equivalent two-mode state 
and evaluating its entanglement. This scaling procedure can 
be iterated to determine all the multimode entanglements ex- 
isting between each mode and each fully symmetric iiT-mode 
sub-block cr^A- (K = 1, • ■ • , N). The 1 x K entanglement 
between the single mode a and any fully symmetric ii'-mode 
partition a-fjK of cr^jv can be determined in a similar way. A 
crucial feature of this scaling structure of the multipartite en- 
tanglement is that, at every step of the cascade, the 1 x K 
entanglement is always equivalent to a 1 x 1 entanglement, so 
that the quantum correlations between the different partitions 
of cr can be directly compared to each other: it is thus possible 
to establish a multimode entanglement hierarchy without any 
problem of ordering. 

To illustrate the scaling structure of multipartite entangle- 
ment in CV systems let us consider a pure, {N + l)-mode 
fully symmetric Gaussian state of the form of Eq. Q. Im- 
posing the constraint of pure state (/i = 1 <=> i^_ = 
i'_i_(jv+i) = 1), one obtains e^ = [1 + fe^(A^ — 1) — 
N - (-l)V(^^ - l)(62(iV + 1)2 - (iV ~ l)2)]/26iV. Such 
a state belongs to the class of CV GHZ-type states discussed 
in Ref. |3]. These multipartite entangled states are the out- 
puts of a sequence of N beam splitters with A^ + 1 single- 
mode squeezed inputs [3, 101 In the limit of infinite squeez- 
ing, these states reduce to the simultaneous eigenstates of the 
relative positions and the total momentum, which define the 
proper GHZ states of CV systems |3|. The CM ct^^hz^ of this 
class of pure states, for a given number of modes, depends 
only on the parameter b = l//i/3 > 1, which is an increas- 
ing function of the single-mode squeezing. Correlations be- 
tween the modes are induced according to the above expres- 
sion for the CO variances e^. Exploiting our previous analysis, 
we can compute the entanglement between a single mode with 
reduced CM /3 and any ii'-mode partition of the remaining 
modes (1 < K < N), by determining the equivalent two- 
mode CM cr'^jiaK- The 1 x K entanglement quantified by the 
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FIG. 1: (color online). Entanglement hierarchy for (A^ + l)-mode 
GHZ-type states (A^ = 9). 



cr'^g.gK- For any nonzero squeezing {i.e. 6 > 1) one has that 

n_ < 1, meaning that the state exhibits genuine multi- 

partite entanglement: each mode is entangled with any other 
iiT-mode block, as first remarked in Ref. |3]. Further, the gen- 
uine multipartite nature of the entanglement can be precisely 



quantified by observing that E' 






> E' 



■m" 



'N 



as shown in 



^m 



logarithmic negativity Ej^ is determined by the smallest 
symplectic eigenvalue fi_ ' 'of the partially transposed CM 



Fig-E The 1x1 entanglement between two modes is weaker 
than the 1x2 one between a mode and other two modes, 
which is in turn weaker than the \ x K one, and so on with 
increasing K in this typical cascade structure. From an op- 
erational point of view, the signature of genuine multipartite 
entanglement is revealed by the fact that performing e.g. a lo- 
cal measurement on a single mode will affect all the other N 
modes. This means that the quantum correlations contained 
in the state with CM cr'^™i can be fully recovered only when 
considering the 1 x TV partition. In particular, the pure-state 
\ X N logarithmic negativity is, as expected, independent of 
N , being a simple monotonic function of the entropy of en- 
tanglement Ev (defined as the von Neumann entropy of the 
reduced single-mode state with CM /3). It is worth noting that, 
in the limit of infinite squeezing (6 — > 00), only the 1 x A^ en- 
tanglement diverges while all the other Ix K quantum corre- 
lations remain finite (see Fig.^. Namely, Ej^ (cr 
-(1/2) log [1 - AK/{N{K + 1) - K{K - 3))], 
not exceed log ^/Z ~ 0.8 for any N and for any K < N. Ai 
fixed squeezing, the scaling with N of the 1 x (TV — 1) entan- 
glement compared to the 1x1 entanglement is shown in Fig.|2l 
(we recall that the 1 x iV entanglement is independent on A^). 
Notice how, with increasing number of modes, the multipar- 
tite entanglement increases to the detriment of the two-mode 
one which becomes distributed between all the modes. We 
remark that this scaling occurs in any Gaussian states, either 
fully or partially symmetric, pure or mixed. For instance, this 
is the case for a single-mode squeezed state coupled with a 
A^-mode symmetric thermal squeezed state. The simplest ex- 
ample of a mixed state with genuine multipartite entanglement 
is obtained from cr'^™i by tracing out some of the modes. 
Fig.|2]can then also be seen as a demonstration of the scaling 
in such a A^-mode mixed state, where the 1 x (A^ — 1) en- 
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FIG. 2: Scaling as a function of iV of tlie 1 x 1 and of the 1 x ( A^ - 1 ) 
entanglement for a (TV + l)-mode GHZ-type CV pure state (b = 
1.5). 



tanglement is the strongest one. Thus, with increasing N, the 
global mixedness can limit but not destroy the genuine mul- 
tipartite entanglement between all the modes. This entangle- 
ment is experimentally accessible by all-optical means |3] and 
it also allows for a reliable (i.e. with fideUty JF > 1/2) quan- 
tum teleportation between any two parties llOll . Therefore, the 
quantification of multipartite entanglement by measurements 
of purity, which, as we have already remarked, can be ex- 
perimentally implemented even without homodyning, leads to 
an accurate estimate of the multi-party teleportation efficiency 
and to direct control on the transfer of quantum information. 

In conclusion, we have shown that multipartite quantum 
correlations of Gaussian states of 1 x A^ partitions under sym- 
metry are endowed with a scaling structure that reduces the 
problem to the analysis of the entanglement of equivalent 
two-mode Gaussian states. Thanks to this reduction, it is 
possible to determine exactly the logarithmic negativity of the 
multimode states and to allow for a reliable experimental esti- 
mate of the multipartite entanglement by direct measurements 
of global and local purities, without the need for the full recon- 
truction of the covariance matrix. Our results apply to many 
cases of practical interest. For instance, the entire class of bi- 
symmetric - i.e. invariant under the exchange of two given 
modes - three-mode Gaussian states ^ has its multipartite 
entanglement completely quantified by the present analysis. 
The generalization of the present approach for the quantifi- 
cation of multipartite CV entanglement to states with weaker 



symmetry constraints and to M x A^-mode partitions (with 
M > 1) awaits further study. 
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Appendix: Proof of the symplectic degeneracy. We prove 
here the multiplicity of the symplectic eigenvalue v^ for the 
CMs cTpN and cr, asserted in Eqs. (|2j and (|9}. We first 
recall that, if E = {vi, . . . ,vm} is the symplectic spec- 
trum of the CM (J, then the 2N eigenvalues of the matrix 
iVtcr are given by the set {^vt}. Let us focus next on the 
CM (Tp'i: in the linear space on which the matrix iVlcrp2 
acts, the eigenvector u_ corresponding to the eigenvalue v^ 
reads u_ = (— i^^^, — 1, z^^^^, 1)^. Due to the symme- 
try of (TpN, any 2A^-dimensional vector v of the form v = 



(0, ...0,-i- 



ei 



4,0. ..0,i- 



ei 



1,0,... 0)"^ (i.e. any 



mode i mode j 

vector obtained by taking w_ in a couple of modes ij and ap- 
pending to it elements for all the other modes) is an eigen- 
vector of iilcTpN with eigenvalue i'_. It is immediate to see 
that one can construct iV — 1 linear independent vectors of the 
above form, proving Eq. (|2j. Clearly, an analogous reasoning 
holds for the matrix cr, proving Eq. (|9jl. 
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